Bianisotropic route to the realization and matching of backward- wave metamaterial 

slabs 
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A concept of backward-wave bianisotropic composite medium matched to free space is suggested. 
It is based on the use of a uniaxial bianisotropic structure embedded into a matrix with negative 
effective permittivity. Since bianisotropy is easier to achieve in the optical range than artificial 
magnetism, this concept is prospective for optical backward-wave metamaterials. As an example of 
possible realizations, a microwave O-composite combined with a wire lattice is analytically studied. 
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Design and studies of materials with negative electro- 
magnetic parameters supporting backward waves is cur- 
rently a very active field of research. The concept of back- 
ward waves is not new: It goes back to the beginning of 
the 19th century and is connected to the names of Lamb, 
Schuster, and Pocklington. In the middle of the 20th cen- 
tury, this concept was extended to waves in homogeneous 
materials and negative refraction effect was theoretically 
predicted (Mandelshtam, Sivukhin, Silin, Veselago). De- 
tailed reviews on the current status and on the history of 
this research field can be found in e.g. [1-7]. 

The main application for backward-wave materials is 
in sub- wavelength imaging devices ("perfect lens" [8]- 
[10]), but a full range of other possibilities is expected, 
especially for optical frequencies, if low-loss and matched 
slabs of backward- wave materials are realized [6, 11]. The 
well-known design approach for the microwave range is 
based on the use of arrays of long thin metal conductors 
and split rings [12]. This approach has been extended 
to terahertz [13] and even infrared frequency range [14]. 
Realization of negative permeability with the use of split 
rings becomes very problematic at optical frequencies, 
and some alternative approaches have been proposed 
[18]-[21]. Considerable progress has been reported along 
this route, but the samples realized so far suffer from 
high losses and poor matching with free space [13, 14]. 
There exist possibilities to realize backward waves also in 
anisotropic media [15, 16] and in anisotropic waveguides 
[17], which do not necessarily require magnetic proper- 
ties of materials, but are restricted to strongly anisotropic 
structures. In addition, difficulties to realize backward- 
wave samples matched to free space inhibit potential ap- 
plications both in the microwave and in the optical re- 
gions. 

Majority of researchers focus on the design of magneto- 
dielectrics, where the backward-wave regime is realized 
when both the permittivity e and permeability fj, have 
negative real parts. Following to paper [22], the bian- 
isotropy is usually considered as a factor that one should 



avoid in the design of backward- wave materials, and ef- 
fort is often concentrated on the design of symmetrical 
variants of split rings which minimizes magnetoelectric 
coupling [23]. However, backward waves can exist in 
more general linear media, namely in bianisotropic media 
(e.g., in chiral media [24]- [26]). It was recently demon- 
strated that in chiral media it is possible to improve the 
characteristics of backward- wave regime [26, 27]. More- 
over, as we will show in this paper, there is a possibility 
to benefit from more design freedom offered by additional 
material parameters. Here we will show that it is possi- 
ble to design a bianisotropic material in such a way that 
it supports linearly-polarized backward waves and a slab 
made of this material is perfectly matched to free space 
for the normal direction of propagation. 

The material can be realized as a composite with 
fi-shaped metal inclusions, and can be called omega- 
medium [28]. Since reliable analytical models of fl- 
particles and fi-media were developed and checked nu- 
merically and experimentally for the microwave range 
[29-31], we use these models (valid below 70 — 100 GHz) 
for a demonstration of the concept having in mind an 
optical realization as one of targets. 

Reciprocal uniaxial f2-media obey the following consti- 
tutive relations [301: 



T> = e-E + j^/e^KJ-H, B=-i-U+j^J^KJ-E 

(1) 

Denoting the unit vector along the optical axis as zq we 
can write the permittivity and permeability dyadics in 
the form: 

e = £o {sJt + £ 2 zoZo) , 71 = Mo UH?t + Mz z oZo) , (2) 

where I t is the two-dimension unit dyadic defined in the 
plane orthogonal to zq: I t = xoXo + yoyo- The magneto- 
electric dyadic is antisymmetric: J = zq x It. Complex 
dimensionless parameter K measures the magnetoelec- 
tric coupling effect. Eigenwaves in such media are lin- 
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early polarized plane waves, similarly to simple magneto- 
dielectrics. 

In this paper we consider a slab of such composite ma- 
terial with the optical axis orthogonal to the slab sur- 
face. For simplicity of the analysis we concentrate on the 
normal-incidence excitation. In this case both eigenwaves 
in the slab have the same propagation constant [30] 



/3 = k y/e t iJ,t - K 2 , 



(3) 



where fco is the free-space wavenumber. For the char- 
acteristic admittances which arc different for the waves 
traveling in the positive and negative directions of the 
axis z (denoted below as Y + and F_, respectively), the 
following relation holds [30] : 



(4) 




where Yq = •v/so/Mo ^ s ^ ne free-space wave admittance. 

In [32] it was shown that omega composites can be 
used to realize absorbing layers which are matched to 
free space. Indeed, one can notice that if the material 
parameters satisfy the condition 



K =2 ^ ~ £ ^> ' 



(5) 



relation (4) gives Y + = Yq. Then the reflection coeffi- 
cient from one side of the slab of an arbitrary thickness 
d equals R = and the transmission coefficient reads 
T = exp(—j(3d). Note that the matching condition can 
be satisfied for arbitrary permittivity and permeability 
values, provided one can control the coupling coefficient 
K. 

If condition (5) is satisfied, the propagation constant 
of the eigenwaves in the medium reads (3) 



= ( £ * • 



(G) 



This result shows that the magnetoelectric coupling not 
only allows one to match the slab to free space, but 
also makes it easier to realize the backward- wave regime. 
From formula (6) we see that the refractive index in a 
matched il-slab is negative when 



Re{e t + fi t } < 0, 



(7) 



which is easier to satisfy than the usual conditions 

Re{e t } < 0, Re{^ t } < 0. (8) 

In particular, it is not necessary to realize negative per- 
meability, which can be more difficult at the optical fre- 
quencies than to satisfy (5). Composites with negative 
permittivity and acceptable losses, on the contrary, are 
available in optics and can be obtained, for example, us- 
ing dilute arrays of metal nanoparticles embedded in a 



dielectric matrix. Like metal split rings [13, 14, 23], metal 
bianisotropic particles can possess a resonance in the in- 
frared range and, possibly, in the visible. Resonant mag- 
netoelectric coupling of these particles can be sufficient 
even if the corresponding resonant magnetic polarizabil- 
ity is not high. This expectation is based on the fact 
that the effect of artificial permeability is of the order 
of /j, ~ O(kod) 2 , while the effect of magneto-electric cou- 
pling is an order of magnitude stronger: K ~ O(kod). 
Here d is the characteristic particle size. 
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FIG. 1: Geometry of the matched backward- wave structure 
and of an individual particle. The uniaxial structure has a 
second identical set of particles lying in the x — z plane (not 
shown), however that second set does not interact with the 
eigenwave with the electric field polarized along y. D c is the 
Q— lattice unit cell size, a m is the period of the wire lattice. 

To study if a realization is possible with a particular 
realistic dimensions we will next consider an example for 
the microwave frequency region. Negative permittivity 
at microwaves can be obtained using a lattice of parallel 
wires (e.g., [33]). A uniaxial omega composite can be re- 
alized as a lattice of pairs of ^-particles. One particle in 
the pair is orthogonal to the other [30, 32]. In this paper 
we restrict our study by the case of the normal propaga- 
tion direction. Then, for the simplified structure shown 
in Fig. 1 all the formulas written above for a uniaxial 
composite are valid. Our purpose is to show that con- 
ditions (5) and (7) correspond to physically achievable 
geometric parameters of f2-particles and wires forming 
the composite medium shown in Fig. 1. 

The microwave electromagnetic model of an omega 
particle represents the particle as a connection of a short 
wire dipole antenna and a small loop antenna [29] . Small 
bianisotropic particles can be characterized by dyadic 
electric and magnetic polarizabilities, which define the 
relations between induced electric and magnetic dipole 
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moments p, m and external electric and magnetic fields 
E, H applied to the particle [28]-[30]: 



p = a ee ■ E + a em ■ H 



m 



E 



H 



(9) 



(10) 



For electrically small particles a lumped-element equiv- 
alent circuit models of such antennas give acceptable ac- 
curacy [30], and the particle polarizabilities can be ex- 
pressed in terms of equivalent circuit parameters. If the 
loop radius a of an fi-particle shown in Fig. 1 is equal to 
or larger than the arm length I, the polarizability com- 
ponents relevant in the case of the normal propagation 
can be approximated as 
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(11) 



(12) 



(13) 



where luq = T = Ri/Lq, and s m is the ma- 

trix relative permittivity. The matrix is assumed to be 
non-magnetic: The matrix permeability \i T = 1. The 
equivalent circuit parameters read 



Ri 



2a ro' 



L = fi a 



log 











Co 



Iog(2//r )' 



(14) 



(15) 



where a is the metal conductivity, and dimensions ro, 
a and / are defined in Fig. 1. Here we have neglected 
the electric polarization of the ring compared to that 
of the dipole portion which is an acceptable approxima- 
tion for electrically small fi-particles [30]. However, the 
bianisotropy of the particle described by formula (13) is 
practically that of the split ring loaded by a lumped ca- 
pacitance Co- This approximation is suitable for the case 
a>l. 

For the relative transversal permittivity and perme- 
ability and for the magnetoelectric coupling the low- 
density approximation (the particles concentration N is 
such that the unit cell of the composite is significantly 
larger than the particles sizes) gives in accordance with 
[32]: 



e t = e r + 



Nail 
£0 



e r + 



so {ujI 



NA 



- UJ 



- j^r) ' 



(16) 



[i t = H r + 



K = Nj 



Na xx 

MO 



fir + 



NBlo 2 



^0 (Wq 



■ juT) ' 



-N 



;D 



{u> 2 - uj 2 + juT) ^/eoMo' 



(17) 



(18) 



where e r and \i r are the relative effective permittivity and 
permeability of the background material, respectively. In 
formulas (16)-(18) we use a different notation for the 
background material permittivity than in (15) because 
in our structure e r ^ e m . The background permittiv- 
ity is, of course, affected by the presence of metal wires 
shown in Fig. 1. However, the capacitance Co is deter- 
mined by the quasi-stationary electric field in the small 
spatial domain of the particle. That is why the matrix 
permittivity e m and not e r enters Co (see more on this 
important issue in [34] and [35]). 

Conditions (5) and (7) should be satisfied at the same 
operational frequency uj. Substituting (16)-(18) into (5) 
we come to equations for the real and imaginary parts of 
(5): 
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(19) 



(20) 



Substituting (16)-(18) into (7) together with (19) and 
(20) it is possible to show that condition e r < — 1 should 
hold for simultaneous validity of conditions (5) and (7). 
This inequality can be satisfied with the negative back- 
ground permittivity created by the wire medium. If the 
wave propagates in the plane orthogonal to the wires and 
the electric field is oriented along them, the wire array 
behaves as a low-frequency plasma. The permittivity of 
the wire lattice is then not spatially dispersive, and we 
have [33, 36] 



£ r (uj) 




(21) 



In [36] one can find two formulas for the plasma frequency 
of such wire arrays, the first one, 



2tt 



P eofioal (log ^- +0.5275)' 



(22) 



is valid if the radius r w of wires is very small (practically 
r w < au,/100). If it is not so small (a/20 < r w < a w /5), 
the second formula is applicable: 



uj 2 = 
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a£,£oMo log 



(23) 
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FIG. 2: Parameters of the bianisotropic layer with the de- 
signed lattice geometry. Effective permittivity (a), perme- 
ability (b), normalized propagation factor (3 /kg (c), reflection 
and transmission coefficients in dB for a 32-mm thick layer 
(d). 



Let us pick up some reasonable values of parameters 
r, s m , a, I, tq (so that a > I and a, I 3> ro) and calculate 
lu . Then from (19) and (20) we find (1 - s r )/N and lo. 
At the frequency lo the medium under design is perfectly 
matched. Then varying N we satisfy the condition e r < 
— 1. Finally we express the needed metal conductivity 
through T: 



a = 



TL r 



(24) 



The results for the particle parameters should satisfy to 
the following physical conditions: 



• Effective medium condition: D c < A e fj/2, where 
A ff = 27r/Re(/3) is the effective wavelength in the 
medium at the working frequency and D c = N~ 
is the size of a unit cell of the O-lattice. 
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• Low concentration condition: D p = max(2a, 21) < 
mm(D c ,a w ). This means that particles do not 
touch each other and the wires. 

• Metal conductivity: a should be that of a known 
metal or alloy. 

If any of these conditions are violated, one can vary the 
initial parameters T, e m , a, I, ro and finally find their 
proper combination. Then we put the wire medium pe- 
riod be equal to the unit cell size a w = D c and find the 
wire medium plasma frequency and the wire radius r w 
from formula (21) together with (22) or (23). These for- 
mulas allow us to study the frequency properties of the 
designed structure. 



A result obtained for the following design parameters: 
e m = 10, I = 3.5 mm, a = 3.9 mm, ro = 0.05 mm, is pre- 
sented in Fig. 2 for the frequency band centered at the 
particle resonant frequency /o = 2.4393 GHz. The per- 
fect matching holds at the frequency / = 2.4359 GHz. 
The parameter T = 0.02 • lo corresponds to the con- 
ductivity a i=a 1.6 • 10 6 1/Ohm-m. With small metal 
losses (those of copper) we could not satisfy the effective 
medium condition. Other parameters were found as ex- 
plained above: N = 2.34 • 10 5 1/m 3 that corresponds to 
— a w rs 16.23 mm and r w — 0.07 mm that allows 
the backward-wave condition over the particle resonance 
band. One can see that the permeability is positive in 
the backward-wave region. The backward-wave regime 
is due to the bianisotropy. The layer thickness was cho- 
sen d = 2D C (the number of particles across the layer is 
enough to use the effective medium approximation) . The 
thickness resonances happen at much higher frequencies 
than / (in the region where e r > 0). 

To conclude, a bianisotropic composite supporting 
backward waves can be obtained with physically real- 
izable parameters of inclusions. This makes the bian- 
isotropic route to matching a backward-wave composite 
be prospective for microwave and hopefully for optical 
metamaterials. Although this matching is perfect for 
the normal incidence only, the mismatch for oblique in- 
cidence should be of the same order as for conventional 
isotropic slabs with [i/e = A*o/ £ o- It is also important to 
notice that the omega slab matching is totally thickness 
independent and holds even for a semi-infinite lattice. 
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